Observation of mixed anisotropy in the critical susceptibility of an ultrathin magnetic 
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Measurements of the magnetic susceptibility of ultrathin Fe/W(110) films with thickness in the 
range of 1.6 to 2.4 ML Fe, show that in addition to the large response along the easy axis associated 
with the Curie transition, there is a much smaller, paramagnetic hard axis response that is not 
consistent with the 2D anisotropic Heisenberg model used to describe homogeneous in-plane ferro- 
magnets with uniaxial anisotropy. The shape, amplitude, and peak temperature of the hard axis 
susceptibility, as well as its dependence upon layer completion close to 2.0 ML, indicate that inho- 
mogeneities in the films create a system of mixed anisotropy. A likely candidate for inhomogeneities 
that are magnetically relevant in the critical region are the closed lines of step edges associated with 
incomplete layers. According to the Harris criterion, the existence of magnetically relevant inhomo- 
geneities may alter the critical properties of the films from those of a 2D Ising model. Experiments 
in the recent literature are discussed in this context. 
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I. INTRODUCTION 

The development of the ability to grow and experimen- 
tally investigate ultrathin ferromagnetic films over the 
last few decades has provided an opportunity to address 
ongoing fundamental questions regarding phase transi- 
tions in two dimensional (2D) systems from a fresh per- 
spective. One such enduring question is the role of mag- 
netic inhomogeneities or defects in altering the behavior 
of the phase transitions of a 2D magnetic system. This 
is important because essentially all physically realizable 
materials have defects. The present article explores these 
issues using magnetic susceptibility measurements of the 
critical region of Fe/W(110) ultrathin films with between 
1.6 and 2.4 Fe ML. 

Ideal ferromagnetic 2D films with in-plane uniaxial 
anisotropy are predicted to display a range of behaviors, 
depending upon the proximity of the temperature to the 
Curie temperature Tq- These films are described by a 
2D anisotropic Heisenberg Hamiltonian, where the crys- 
talline and surface anisotropies are most often modelled 
as an exchange interaction that is anisotropic at the level 
of 1% to 0.1%. Both theory^ and simulations 2 indicate 
that the presence of this small anisotropy moves the Curie 
temperature from K to of order 100s of K, and induces 
the critical behavior of a 2D Ising model. A growing num- 
ber of experimental measurements of the critical expo- 
nent of the magnetization, /3, confirm this^ Simulations 4 
indicate that the cross-over from 2D anisotropic Heisen- 
berg to 2D Ising behavior occurs relatively far from Tc, 
and, to our knowledge, has not yet been observed ex- 
perimentally. The hard axis susceptibility in the criti- 
cal regime can also be derived from the 2D anisotropic 
Heisenberg model. 5 It is a measure of the microscopic 
anisotropy that persists in the paramagnetic region, and 
is described qualitatively by the 2D Ising models 

Real films have structural inhomogeneities that can 
easily create local magnetic inhomogeneities. Examples 



are missing atoms (site dilution), structural defects or 
atomic adsorption that can locally alter the exchange 
coupling or the microscopic anisotropy, or both. The spe- 
cific example of a sample where sites have one of two dif- 
ferent anisotropies is termed "mixed anisotropy" X Cen- 
tral questions are whether these inhomogeneities are rel- 
evant or irrelevant to the critical behavior, and how they 
affect the extent of the asymptotic region. The Harris 
criterion^ treats the case of random, point-like inhomo- 
geneities. Experiments on 3D systems have confirmed 
the Harris criterion.—^ 1 In particular, for 3D Ising sys- 
tems, introducing inhomogeneities in the form of either 
altered exchange coupling along the Ising axis^, or local 
XY anisotropy perpendicular to the Ising axisii gives al- 
tered critical exponents in agreement with the 3D site di- 
luted Ising models In two dimensions, inhomogeneities 
become relatively more important. In the case of the 2D 
Ising model, the Harris criterion is not predictive. This 
has lead to a large body of theoretical and numerical 
studies of the effect of inhomogeneities in this system^ 

If the inhomogeneities are not point-like and random, 
but correlated to form a line, then it may be that the 
effective dimensionality of the system is changed^ so 
that the prediction of the Harris criterion is altered<&i£~— 

The current study investigates Fe/W(110) ultrathin 
films near the Curie transition, using measurements of 
the magnetic susceptibility. This system has very strong 
uniaxial, in-plane anisotropy 1 - and is a model for 2D 
Ising behavior. An experimental study of the magne- 
tization in the ferromagnetic phase by Back et al^ has 
shown 2D Ising power law scaling over 18 orders of mag- 
nitude in reduced variables. In addition to the divergent, 
easy axis response, the present work reports a small re- 
sponse in the hard axis, in-plane magnetic susceptibil- 
ity in the paramagnetic phase. The form of the hard 
axis response and its systematic dependence on the film 
thickness and layer completion demonstrates that it is 
due to inhomogeneities that create a system with mixed 
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anisotropy. A likely candidate for the relevant inho- 
mogeneities is the closed lines of atomic steps that are 
present when the film thickness deviates from two com- 
plete monolayers. The implications of these magnetic in- 
homogeneities for the critical behavior is discussed with 
reference to recent experimental findings of the variation 
of the critical exponent of the susceptibility, 7, with layer 
completion in Fe/W(110)3 



II. EXPERIMENTAL RESULTS 

Ultrathin Fe films grown on W(110) in ultrahigh vac- 
uum have been studied intensively, and have revealed 
magnetic behavior that is sensitive to film thickness, 
structure and cleanliness. The first monolayer of iron 
wets the W(110) surface very well^i and exhibits an in- 
plane crystalline anisotropy along [1-10] The interfa- 
cial anisotropy of 1.9 mJ/m 2 is one of the largest known 
for an ultrathin film system. According to scanning tun- 
neling microscopy (STM) images, the structure of the 
second monolayer depends upon the substrate orienta- 
tion and temperature during growth. Deposition at or 
near room temperature on a substrate cut at an angle 
< 0.3° to the [110] surface normal produces islands of 
2ML thicknes s 19 ! 22 " — which grow and begin to coalesce 
near a thickness of 1.6 ML. There is good layer com- 
pletion near 2.0 ML, with small regions of 1 ML and 3 
ML Fe coexisting. If the substrate is miscut at an angle 
> 0.5°, and the Fe is deposited near 660 to 700K, then 
the presence of many atomic steps and the increased mo- 
bility of the Fe atoms induces step-flow growth ) 21 i 25 i 26 
resulting in long "strips" of 2 ML thickness parallel to 
the substrate steps, alternating with strips of 1 ML or 3 
ML. 

In the present experiments, the crystal miscut is < 0.2° 
as determined using ex situ STM imagesi^l The crystal 
was cleaned according to standard procedures of alter- 
nate heating in 1 x 10~ 7 Torr oxygen, and flashing to 
white heat until the surface carbon contamination was 
not detectable using Auger electron spectroscopy, and a 
sharp W(110) LEED pattern was obtained. The first 
monolayer was deposited at room temperature and an- 
nealed for 150 s at 500 K. A second deposition of between 
0.6 and 1.4 ML to complete the film was made at room 
temperature. LEED images confirmed that the Fe films 
were pseudomorphic with the substrate in this range of 
coverages 2 ^. Although in situ STM imaging was not pos- 
sible, layer growth through the coalescence of monolayer 
islands is expected based upon the previously cited stud- 
ies. In order to confirm this expectation, the film growth 
was monitored using Auger electron spectroscopy and 
pseudomorphic growth was confirmed using low energy 
electron diffraction. Fig. la presents the attentuation of 
the W(110) Auger signal as Fe was evaporated onto it in 
a series of sequential depositions. After each deposition, 
the film was annealled to 700K. The clear break in the 
plot illustrates the completion of one monolayer, where 
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FIG. 1: a) The film thickness was calibrated by plotting the 
attenuation of the W(UO) substrate signal of a Fe film grown 
in sequential steps against the deposition time. The film was 
annealed at each step. The break in the curve indicates a 
thickness of 1 ML and good wetting of the substrate, b) The 
W Auger attenuation is plotted against Curie temperature 
(peak in X[i-iol) f° r the films in this study. The literature 
value of the Curie temperature for 2 ML Fe/W(110) corre- 
sponds to an attenuation of 0.38±0.01, in very good agree- 
ment with the expected value for layer-by-layer growth of 2 
ML, 0.36±0.024. 



the W(110) Auger signal is attenuated to 0.60 ± 0.02 of 
the value for a clean substrate. 

Fig. lb plots the Auger attenuation coefficient of 
the W substrate upon deposition of the entire Fe film, 
as a function of the Curie temperature, for the films 
used in this study. Layer growth of a second complete 
monolayer is expected to correspond to an attenuation 
of (0.60 ± 0.02) 2 = 0.36 ± 0.024. According to the 
literatures^ the Curie temperature of 2 ML Fe/W(110) 
is 455±3K. In fig. lb, the W Auger signal is in fact atten- 
uated by 0.38±0.01 at this temperature, independantly 
confirming very good layer completion at 2 ML. The fact 
that the W signal is not attenutated quite as much as ex- 
pected indicates a residual small net area of 1 and 3 ML 
regions. The magnetic properties, as indicated by Tq, are 
consistent with previous studies where island growth has 
been demonstrated by STM^ Since the Curie tempera- 
ture is a much more sensitive measure of relative thick- 
ness than the Auger attenuation, Fe thicknesses quoted 
in the remainder of this article are calculated based on 
the Curie temperature and the calibration curves in fig. 
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FIG. 2: Real (solid dots) and imaginary (open dots) 
parts of the magnetic susceptibility of 1.7 ML Fe/W(110). 
a) X[ i-io](T) b) X [ooi](T) 



1. 

Previous studies^ have shown that isolated 2 ML Fe 
islands on 1 ML Fe/W(110) have perpendicular mag- 
netic anisotropy at temperatures below about 235 K, but 
that once the islands coalesce at a thickness of about 
1.6 ML, there is once again strong in-plane, uniaxial 
anisotropy. At lower temperatures, the perpendicular 
anisotropy persists to higher coveragei 24 i 30 The perpen- 
dicular anisotropy of the islands is very sensitive to gas 
adsorption^ and a small gas exposure induces a spin- 
reorientation transition that leaves the entire sample with 
in-plane anisotropy. Since the films in the present study 
were all of thicknesses greater than 1.6 ML Fe, and the 
temperature range investigated was 400 - 480 K, only 
in-plane anisotropy was observed, consistent with these 
results. 

Measurements of the magnetic susceptibility were 
made in situ using the surface magneto-optic Kerr ef- 
fect with an applied ac field of 150 Hz and lock- in 
amplification4i Recent modifications of the apparatus 32 
that increased structural rigidity, reduced scattered light 
from the polarizing crystals, and dispensed with aper- 
tures, have increased the sensitivity and reduced the 
noise so that optical rotations of 15 nrad are detectable. 
The measurements use a pair of field coils aligned very 
nearly parallel to the sample surface, and in the scatter- 
ing plane of the laser light. The light is incident at 45° 
to the sample surface, so that the measured ellipticity 
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FIG. 3: a) The FWHM of the susceptibility peak is plotted 
as a function of the applied ac field amplitude, for three ori- 
entations of the scattering plane and field direction; [1-10], 
[001] and between these orientations, 10° from the [001]. b) 
The maximum amplitude of the real and imaginary parts of 
the susceptibility are plotted as a function of the angle that 
both the scattering plane and direction of applied field make 
with the [001] hard axis. 



is a combination of the longitudinal and polar Kerr ef- 
fects. The sample can be rotated about its normal, so 
that the aligned scattering plane and applied field can 
lie along any in-plane direction. For definiteness, the 
measured susceptibilities are labelled with a subscript 
that denotes this in-plane direction. The susceptibility 
is given in units of /iradians/Oe, since these can be cali- 
brated absolutely. 

Figure 2 presents measurements of the magnetic sus- 
ceptibility for a 1.7 ML Fe/W(110) film using an ac field 
of 0.50 Oe. Part a) contains the real and imaginary parts 
°f X[i-io](T)- The narrow peaks in both the real and 
imaginary easy axis response are characteristic of a Curie 
transition and the associated hysteretic dissipation just 
below Tq ■ This result is entirely as expected for remanent 
magnetization along the easy axis. Part b) illustrates the 
hard axis susceptibility XrooilC^)- This curve is different 
from the easy axis result in three important ways: the 
amplitude is smaller by a factor of 15, the full width at 
half maximum (FWHM) of the peak is significantly nar- 
rower, and there is no measurable imaginary response. 

A systematic study of this film as a function of the crys- 
tal orientation is presented in figure 3. Part a) presents 
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the FWHM of the real part of the susceptibility as a 
function of the applied ac field measured along the easy 
[1-10] axis, the hard [001] axis, and an intermediate in- 
plane angle 10° from the hard axis. These results make 
it clear that the peak widths are not due to finite field ef- 
fects, and quantifies the fact that the hard axis response 
is narrower than the easy axis response by nearly a factor 
of 2. Fig. 3b presents the amplitude of the real and imag- 
inary parts of the susceptibility as a function of the angle 
<t> that the in-plane hard axis makes with the scattering 
plane and applied field. Both of these amplitudes are well 
described by a quadratic variation about the hard axis, 
where the imaginary response goes to zero. This confirms 
that the measured hard axis response is a distinct signal 
that has no imaginary part. Moving away from the hard 
axis creates a component of the applied field along the 
easy axis, and the measured susceptibility is a mixture 
of the easy and hard axis responses that is dominated by 
the much larger easy axis susceptibility. 

All of the above experiments were performed with p- 
polarized light. Measurements with s-polarized light (not 
shown) produced similar susceptibility traces. This ex- 
cludes the possibility that the signals labeled X[ooi] arise 
through the transverse Kerr effect with the magnetiza- 
tion perpendicular to the scattering planed Measure- 
ments were also made with a separate magnetic coil 
aligned very nearly with the surface normal hard axis, 
[110] (not shown). These signals were a further order 
of magnitude smaller than those along the in-plane hard 
axis, and had the shape of the easy axis signal when the 
scattering plane included the [1-10] and [110] direction, 
and the shape of the in-plane hard axis signal when the 
scattering plane included the [001] and [110] directions. 32 
It was therefore concluded that these are not true mea- 
surements of the susceptibility normal to the surface, but 
rather small components of the in-plane susceptibilities 
due to a small misalignment of the coil and crystal axes 
that could not be removed using the degrees of freedom 
of the crystal holder. 

The results of fig. 2 and 3 indicate that these films 
have a reproducible, narrow hard axis response that con- 
tains no imaginary component - i.e. it is dissipationless. 
In order to better characterize the system, susceptibility 
measurements have been made on a collection of 22 sep- 
arate films with thicknesses in the range of 2.0 ±0.4 ML 
of Fe. The main panel of figure 4a presents three curves 
measured sequentially from a 2.2 ML film. First, the 
real part of the easy axis susceptibility (solid circles) was 
measured, then the sample was rotated and the hard axis 
susceptibility (solid triangles, and in more detail in the 
inset) was measured. The sample was then rotated back 
to its original alignment and the easy axis susceptibility 
was remeasured (open circles). For this film the peak 
of the hard axis susceptibility occurs about 3K higher in 
temperature than the peak of the easy axis susceptibility. 
This is not because of thermally induced changes in the 
film, since the easy axis signal is reproducible to within 1 
K after measuring the hard axis curve. This reproducibil- 
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FIG. 4: a) The real part of the susceptibility of a 2.2 Fe ML 
film measured three times in sequence: along the easy axis 
(solid circles), the hard axis (solid triangles, and inset), and 
along the easy axis again (open circles). The hard axis signal 
peaks about 3 K above the easy axis signal. The latter is re- 
produced to within 1 K. b) The hard axis [001] susceptibility 
is plotted for films of different Fe thickness, with a tempera- 
ture scale measured relative to the peak temperature of the 
corresponding easy axis [1-10] susceptibility. The uncertainty 
in this relative scale is ± IK. 



ity was checked on 1/3 of the films, either by sequential 
easy-hard-easy axis measurements or by repeated mea- 
surements of the easy or hard axis susceptibilities. This 
includes samples showing hard axis susceptibility peaks 
up to 10 K above Tc- An uncertainty of ±1K of the 
peak positions is a reliable representation of our mea- 
surements, as no instance of a larger change in Tc upon 
re-measurment was observed. 

Fig. 4b shows hard axis [001] susceptibility of a se- 
lection of films, labelled by their thickness. Note that 
the temperature scale is relative to the peak of the easy 
axis susceptibility for each film, T e , which closely tracks 
Tc- The uncertainty in this relative scale is ±1 K, as 
discussed above. These curves illustrate that the hard 
axis response depends upon the film thickness in two im- 
portant ways: first, the amplitude of the susceptibility 
decreases as the film thickness moves closer to 2.0 ML; 
second, there is a correlated shift in the peak tempera- 
ture, so that the hard axis peak moves further into the 
paramagnetic phase as the film approaches a thickness of 
2.0 ML. In all cases, there is no imaginary response. 

Figure 5 summarizes the results for the entire collec- 
tion of films. Part a) plots the amplitude of the hard 



5 



Fe thickness (ML) 

1.66 1.87 2.08 2.30 2.51 




O 420 440 460 480 500 

Q_ 

Peak temperature T e (K) 

FIG. 5: Features of the hard axis susceptbility for 25 differ- 
ent films are summarized, a) The peak amplitude is plotted 
against T e . The thickness scale on the top of the plot is de- 
rived using the data in fig. 1. b) The peak position, relative 
to the peak temperature of the easy axis susceptibility of the 
same film. The non-linear right hand scale uses the reduced 
temperatue t = Th/T e — 1. The fitted curves are discussed in 
the text. 



axis susceptibility as a function of T e . Part b) plots the 
relative peak position of the hard axis response, Th, com- 
pared to that of the easy axis susceptibility, T e . The error 
bars indicate the characteristic reproducibility upon re- 
peated measurement, as was discussed previously. These 
data confirm that the hard axis susceptibility is part of 
the paramagnetic, non-hysteretic response of the system 
as the incipient ferromagnetic state develops with the ap- 
proach to Tq from above. Both the amplitude and the 
peak temperature of the paramagnetic hard axis response 
are strongly correlated with Fe coverage, 9, as some func- 
tion of 1 6 - d 1, where 6q= 2.0 ML. The fitted lines are 
discussed in the next section. The symmetry of the data 
about a thickness very close to 2 ML indicates that re- 
gions of 1 ML and 3 ML in a nominally 2 ML film produce 
equivalent effects in the hard axis susceptibility. 



III. DISCUSSION 

The form of the hard axis susceptibility and its de- 
pendence upon \9 — 8q I indicates that these films are 
not homogeneous, anisotropic 2D Heisenberg systems. 
Simulations*^ theory and experiment^ indicate that in 



this model the hard axis susceptibility exhibits a very 
small, rounded, non-divergent contribution that is maxi- 
mized at Tc and persists below the Curie temperature, in 
qualitative agreement with the 2D Ising model. 5 Strongly 
anisotropic behavior persists relatively far from Tc~ In 
contrast, the present experiments observe a hard axis sus- 
ceptibility that is sharply peaked in the paramagnetic re- 
gion above Tc, does not persist down to Tc, and is only 
1 to 2 orders of magnitude smaller than the easy axis 
signal. In addition, the susceptibility of a system that is 
effectively uniform should not vary as a strong function 
of|0-0 o |. 

The measurements are consistent instead with a sys- 
tem of mixed anisotropy. Magnetic inhomogeneities that 
can produce mixed anisotropy arise naturally from the 
structural inhomogeneities implicit in films that are not 
complete monolayers. In a system such as Fe/W(110), 
where layer-by-layer growth near 2 ML proceeds by the 
colescence of islands, the distribution of 1 ML pits and 
3 ML islands will vary in number, perimeter and area as 
some function of \9 — 9q\. Changes in the magnitude 
and/or direction of anisotropy due to monolayer step 
edgesj^l atomic adsorption at step edges^ or a mono- 
layer change in thickness^ are well known. All of these 
effects can give rise to a minority of the sample which 
has a local easy axis which is not aligned with the global 
easy axis of the 2D anisotropic Heisenberg system. The 
local easy axis response from these regions would then 
appear in the global [001] hard axis susceptibility. 

The functional form of the dependence of the paramag- 
netic hard axis susceptibility in fig. 5 can be understood 
qualitatively in terms of inhomogeneities of characteris- 
tic size d, by considering the changing correlation length, 
£, of the system as the temperature approaches Tc- At 
high temperatures, the correlation length will be small, 
such that £ << d. In this case there will be essentially 
no coherent magnetic response from different regions of 
the inhomogeneity, and the hard axis susceptibility will 
be near zero. At some temperature closer to Tc, the 
condition £ « d will be met. The entire inhomogeneity 
can then respond coherently and produce a maximum in 
the hard axis susceptibility. Even closer to Tc, £ >> d. 
In this circumstance, the inhomogeneity is immersed in 
a large correlated region with a uniaxial easy axis. The 
resulting exchange stiffness within the correlated region 
inhibits the hard axis response of the defect atoms, and 
the hard axis susceptibility falls quickly once more. This 
creates a peak in the hard axis susccptibililty. 

According to this qualitative picture, the temperature 
of the peak of the hard axis susceptibility, Th , will occur 
when £ w d. Since the coverage 6 represents a surface 
area on the film, the characteristic size or diameter of 
the inhomogeneity will scale as d ~ W\6 — 6q\, giving 

Th - T e ^ T h - T c , Co y I — (i) 

T e ~ T c ^ { Z{T h ) } ~ d ~ y/ie-Gol' 

since v = 1 for the 2D Ising model. At temperature Th, 
the maximum value of the hard axis susceptibility will 
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scale as the total number of atoms in a characteristic 
inhomogeneity. This leads to scaling as d n , where the 
(integer) n depends upon the geometry of the inhomo- 
geneity. 



x(T h ) 



(2) 



In these expressions, A and B are constants. For the 
fitted lines in fig. 5, the primary, more precise variable 
\J\T e — T°| has been used instead of the derived vari- 
able yJ\e -6 \. The fitted values for T e ° are 453.0±0.7 
K and 454.7±0.5 K, in fig. 5a and 5b respectively. Both 
values are very close to the Curie temperature of 2 ML 
Fe/W(110). The best fit to the data in fig. 5a is given 
with n = 1. 

It is not possible to unambiguously identify the inho- 
mogeneities responsible for the hard axis susceptibility 
in the present experiments, but the most likely candi- 
date is the monolayer steps at the edges of regions of 1 
or 3 ML Fe thickness. Near a coverage of 2 ML, these 
step edges will form closed curves along the perimeter of 
3 ML islands or 1 ML pits of characteristic size d that 
will contribute to the susceptibility as ~ d 1 . This agrees 
with the observation of identical magnetic effects for a 
local monolayer increase or decrease in thickness from 2 
ML, since the creation of an island or pit gives an equiv- 
alent perimeter. A less likely candidate is the areas of 
the 3 ML islands and I ML pits. These do not normally 
have equivalent magnetic effects, as can be seen, for ex- 
ample, in the monotonic change in the Curie temperature 
with thickness. In addition, areas scale as n = 2 in eq.©. 
This gives a linear fit to a small range of data close to 2.0 
ML, but does not represent all the data nearly as well. It 
is unlikely that the observed changes in the anisotropy 
are due to adsorption, since then they would not de- 
pend upon layer completion, as observed, but rather the 
elapsed time. 

Although it is obvious that there must be disorder at 
some level in a real system, these results demonstrate 
that the disorder related to layer completion is magneti- 
cally relevant and creates measurable, quantifiable effects 
in the critical region. Since the inhomogeneities create 
an additional response along a different axis than the 
bulk of the film, the disorder creates a film with mixed 
anisotropy. This highlights the sensitivity of the para- 
magnetic susceptibility to inhomogeneities, since these 
measurements permit an observation of their explicit 
magnetic response. 



IV. IMPLICATIONS FOR CRITICAL 
BEHAVIOR 

These results are strong evidence that the Fe/W(110) 
films have mixed anisotropy, with the degree of disorder 
depending sensitively on the degree of layer completion 
near 2 ML. This being the case, there are very general ar- 
guments concerning how these magnetic inhomogeneities 



will affect the critical properties of the films, that apply 
even though the nature of the magnetic inhomogeneities 
is not known precisely. There are three inter-related as- 
pects to consider. The first is the local change in the 
magnitude of the anisotropy at an inhomogeneity and the 
resulting local variation in the Curie temperature. The 
second is the characteristic size of the inhomogeneities, 
and the third is the overall degree of disorder created by 
the inhomogeneities. 

Associated with the distribution of anisotropy direc- 
tions and magnitudes produced by the inhomogeneities, 
there will also be a distribution of local Curie temper- 
atures. This is because the the Curie temperature of 
the anisotropic 2D Heisenberg model is determined by 
the anisotropy, and would be zero in the absence of an 
anisotropyJ^ An estimate of the Curie temperature for 
an anisotropic 2D Heisenberg system with exchange cou- 
pling J and anisotropy K is given by Serena et 
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In this expression, Tq is the Curie temperature of the 2D 
Ising model with exchange coupling J, and Tq is the 
Curie temperature predicted by a first order spin wave 
treatment of the anisotropic 2D Heisenberg model. The 
latter is given by 



k B T™ = 



4ttJ 

MjJk)' 



(4) 



For a small change in the anisotropy, the change in the 
Curie temperature is 



AT, 



c 



AK 



T c (2tt/2.27) + ln( J/K) K 



(5) 



In the current system with an anisotropy of Kj J w 0.001 
and Tc = 450X, a change in the Curie temperature 
of about 5K can be produced by a 10% change in the 
anisotropy. 

The Harris criterion^ considers the question of whether 
or not these local variations in the Curie temperature will 
change the critical properties, for instance the critical ex- 
ponents, of a system with inhomogeneities. If the inho- 
mogeneities are not point-like, but correlated, then the 
Harris criterion finds them relevant and able to change 
the critical properties ^ 15 ' 17 The characteristic size of the 
inhomogeneity is important, since, in practical terms, 
an inhomogeneity of size d is point- like if £ >> d. If 
£ < d, then the inhomogeneity represents an indepen- 
dent region that contributes to the measured suscepti- 
bility. The present experiments observe this independent 
contribution in the hard axis [001] susceptibility, where it 
is prominent because the contribution from regions with 
[1-10] easy axis anisotropy is so small. The existence 
of this signal is an indicator that £ < d, that the homo- 
geneities are not point-like in this temperature range, and 
that critical exponents determined in this range are not 
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likely to reflect universal properties. The scaling of the 
peak temperature and amplitude with y/9 — 6q indicates 
an increase in d with layer disorder and a reduced range 
close to Tq where the inhomogeneities can be considered 
point-like. 

Finally, even if the inhomogeneities are point-like and 
randomly distributed, the overall degree of disorder can 
still be important. The Harris criterion states that the 
relevance of inhomogeneities depends upon the sign of 
the critical exponent of the specific heat, a. Since the 2D 
Ising model has a = 0, the Harris criterion is not predic- 
tive. However, theoretical 3 ^— and computationa l 13 ' 39 
work finds that the critical exponents are unchanged, 
but that additional logarithmic terms become important 
when the reduced temperature t = T/Tq — 1 < t cr , where 
lat cr = — 1/g depends upon the overall degree of disor- 
der through the parameter g. The parameter g — > as 
disorder disappears, but is not otherwise absolutely quan- 
tified except through numerical studies. Numerical work 
by Roder et al.« 39 has indicated that fitting to a power law 
when logarithmic corrections are important produces 



7e// ~ HsingW + 



1 ln(|lni| 
2 
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In this case j e f f is always greater than the 2D Ising value. 

In this context, the present experimental results bring 
new insight to a recent study-^ of the critical exponent, 
7, of the magnetic susceptibility of Fe/W(110). In this 
study of 25 films, in-plane easy axis ac magnetic suscep- 
tibility measurements have been fit objectively to deter- 
mine four critical parameters simultaneously by the min- 
imization of the variance from a simple power law in In \ 
vs. Int. In addition to the intercept xo an d slope 7 e //, 
both Tq and a fitting cutoff In t x have been found. t x de- 
fines the minimum of the range of reduced temperature 
over which the ln-ln plot is linear. It defines how close 
to Tc the power law extends before either extrinsic or 
intrinsic effects cause the susceptibility to deviate from a 
power law. The maximum range of reduced temperature 
that is available, t ma x, is determined by experimental 
signal-to-noise considerations. 

The data from ref.[20] is replotted in two ways in fig. 
6. Part a) plots the fitted value of j e ff against the peak 
temperature, T e , of the easy axis susceptibility (essen- 
tially the Curie temperature). These data show that the 
2D Ising value of 7 is recovered consistently and repro- 
ducibly only when T e ss 455K, which corresponds to the 
Curie temperature of 2.0 ML Fe films. In fig. 6b, the 
same data is plotted against the cut-off lni^. \nt max is 
indicated by the dashed line. Points further to the left on 
the plot represent susceptibility measurements that dis- 
play power law scaling extending closer to Tc ■ This plot 
shows two interesting correlations. First, the Fe films 
where Jeff takes the 2D Ising value are not only very 
near 2 ML in thickness, they are also films which have 
the longest range of power law scaling. Second, for films 
where 7 e // is not the 2D Ising value, there is a systematic 
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FIG. 6: A summary of experimental results from refi2p| for 
j e f f , determined from measurements of the in-plane easy axis 
magnetic susceptibility of Fe/W(110) films, a) The value of 
the exponent is plotted against the peak temperature, T e ~ 
Tc- 2D Ising values are obtained in a narrow temperature 
range corresponding to a Fe thickness of very nearly 2 ML. b) 
The same data is plotted against lnt x , where t x is the fitted 
minimum limit of power law behavior on a ln-ln plot. The 
maximum limit of power law behavior, due to signal levels, is 
indicated by the dashed line at \nt ma x- 



relationship such that the value of j e f / increases as the 
range of power law scaling is reduced. 

For Fe/W(110) films that are very nearly 2 ML in 
thickness, comparison of fig. 6 and fig. 5b (using the right 
hand scale in In t) indicates that the inhomogeneities give 
no signal in the hard axis susceptibility within the fit- 
ting region for Jeff, and are therefore point-like in this 
range. Nonetheless, the peak in the hard axis suscepti- 
bility at higher temperature shows that there is disorder, 
and the parameter g is non-zero. The fact that the fitted 
value of 7e/ / for the 2 ML films is that of the 2D Ising 
model indicates that t cr < t x , and the logarithmic cor- 
rections do not affect the fitting. This establishes a limit 
of g < -l/m(i x ) = 0.15 for the 2 ML films. This is mid- 
way between what numerical studies calibrate as strong 
disorder-^ (g = 0.30) and weak disorder^ (g = 0.017). 

For films that depart even 0.1 - 0.2 ML from a complete 
2 ML, the situation is more complicated. The peak in 
the hard axis susceptibility in fig. 5 moves to hit = —4.5, 
which is midway through the fitting region for j e f f in fig- 
6. In this circumstance, the inhomogeneities may not be 
point-like, and variations in the local Curie temperature 
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may not be averaged out. At the same time, the increased 
disorder will increase g, moving t cr to the right. If In t cr > 
lnix, the eq.© indicates that logarithimic corrections 
cause Jeff ~ 2.0. Since there is no independent method 
of determining g for the films, this possibility cannot be 
tested. By some combination of these two effects, the 
fitted value j e f / departs from the result for the 2D Ising 
model. 

As the film thickness departs even further from 2 
ML, the peak temperature of the hard axis susceptibil- 
ity moves close to t x , so that in the entire fitting region 
the inhomogeneities are not point-like. The temperature 
range where they are point-like moves so close to Tq that 
it is inaccessible to the experiment, and the question of 
logarithmic corrections does not arise. 

The previous discussion indicates how the character- 
istic size of the defects causes the experimentally deter- 
mined critical exponent to deviate from the value of the 
2D Ising model, but it does not indicate why the devia- 
tion is always to a larger value, and why it is correlated 
to lnt x in fig. 6b, but not \6 — 9 \ in fig. 6a. In ref . (|20T ) . 
the authors state that a systematic increase in Jeff as 
seen in the experiment is mimicked by generating "data" 
numerically using a 2D Ising model with a Gaussian dis- 
tribution of Tc of width 2 to 5 K. As the range of Tc is in- 
creased, so does 7e/ / • This is entirely consistent with the 
present findings, where inhomogeneities create a mixed 
anisotropy. According to cq.©, a modest 5% to 10% 
change in the local anisotropy at the inhomogeneity can 
create the required range of Curie temperatures. If the 
characteristic size of the inhomogeneities is small enough 
(very close to 2 ML), then the inhomogeneities are point- 
like in the critical region and the range of local values of 
Tc is not relevant. However, moving away from 2 ML, 
the inhomogeneities are not point-like in the temperature 
region that is fitted for j e f f , and the range of anisotropy 
values determines the range of Tc and affects the fitted 
value of Jeff- Independent films of the same thickness 
can have a different distribution of anisotropies that pro- 
duce a range of precisely determined, but different values 
of Jeff, as is seen in fig. 6a. However, regardless of thick- 
ness, the larger the range of Tc in each film, the shorter 
will be the range in reduced temperature that will appear 
linear on a In- In scale. Therefore, \nt x is correlated with 
jeff as in fig. 6b. 

There are very few other published studies with which 
to compare these results. Most critical studies of ultra- 
thin magnetic films measure the exponent of the magne- 
tization, /?, in the ferromagnetic state below the Curie 
temperature^ In this case, the ferromagnetic correlation 
length is very large and the inhomogeneities are always 
point-like. We are unaware of any other systematic stud- 
ies of the critical susceptibility as a function of layer 
completion. However, two studies of individual films are 
noteworthy. Jensen et aL- investigate a Co film on a vic- 
inal Cu(l 1 17) surface, and find only a small hard axis 
susceptibility consistent with a homogeneous anisotropic 
2D Heisenberg model. The regularly-spaced lines of step 



edges of vicinal surfaces provide the dominant anisotropy, 
and do not form closed curves of characteristic size d. 
We are unaware of measurements of critical exponents of 
films grown on this type of surface, but if theoretical cal- 
culations that consider the effect of infinite line defects 
in the 2D Ising model are applicable ; 14 ' 15 ' 17 they may be 
nonuniversal. 

Finally, the results of Back et aZ. 19 seem to contradict 
our findings. They study a 1.7 ML film where the first 
monolayer is not annealed and the resulting Tc ~ 340K. 
However, the film has an island structure, and, in a sep- 
arate publication 4 ^, they find a broad in-plane hard axis 
susceptibility centered near Tc when a relatively large 
field of 12 Oe is applied. The fact that they measure the 
2D Ising value of j may be a result of the use of a static 
method of measurement that can more closely approach 
Tc than our ac method^ 



V. CONCLUSIONS 

Measurements of the magnetic susceptibility of 
Fe/W(110) films in the range of 1.6 to 2.4 ML Fe have 
revealed a small, hard axis response that is distinct from 
the much larger easy axis susceptibility associated with 
the phase transition to ferromagnetism. The shape, size 
and peak temperature of the hard axis susceptibililty is 
not compatible with the 2D anisotropic Heisenberg model 
that is generally considered to describe a homogeneous, 
uniaxial, in-plane ferromagnet. In addition, the depen- 
dence of the signal on the layer completion of the film 
as \9 — 8 Q \, where 9 — 2.0 ML, is not consistent with a 
homogeneous system. Rather, the hard axis signal in- 
dicates that magnetic inhomogeneties that remain rele- 
vant in the critical region create a system with mixed 
anisotropy. While the majority of the system can be de- 
scribed by the 2D anisotropic Heisenberg model (2D Ising 
in the critical region), a minority display a local easy axis 
behavior along the global hard axis. Although the pre- 
cise nature of the inhomogeneities has not be identified, 
a likely candidate is the closed lines of monolayer step 
edges at the perimeter of the islands or pits associated 
with an incomplete layer. 

Because the system has mixed anisotropy, the mag- 
netic inhomogeneities affect the critical behavior of the 
film. The variation of the peak amplitude and temper- 
ature of the hard axis susceptibility with \/8 — 9q indi- 
cates that the inhomogeneities have a characterisitic size 
d. Since the peak occurs when £ tzs d, it marks a tempera- 
ture range where the inhomogeneities are not point-like. 
In this case, the local variations in Tc induced by the 
local changes in anisotropy are not averaged over, and 
lead to an increase in the value Jeff fitted to the sus- 
ceptibility in a previous study. Only when the Fe films 
are very close to a complete 2 ML is d small enough that 
the inhomogeneities are point-like over the temperature 
range accessible to the experiments, and the fitted expo- 
nent has the 2D Ising value. Although this makes it more 
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difficult to access the asymptotic region experimentally, 
it offers an interesting opportunity to study the explicit 
magnetic response of defects and inhomogeneities. 
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